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Abstract 

Pseudoheating of ions in the presence of Alfven waves is studied. We show that this 
process can be explained byExB drift. The analytic solution obtained in this paper 
are quantitatively in accordance with previous results. Our simulation results show 
that the Maxwellian distribution is broadened during the pseudoheating; however, the 
shape of the broadening distribution function depends on the number of wave modes 
(i.e., a wave spectrum or a monochromatic dispersionless wave) and the initial thermal 
speed of ions (v p ). It is of particular interests to find that the Maxwellian shape is more 
likely to maintain during the pseudoheating under a wave spectrum compared with a 
monochromatic wave. It significantly improves our understanding of heating processes 
in interplanetary space where Alfvenic turbulences exist pervasively. Compared with a 
monochromatic Alfven wave, E x B drift produces more energetic particles in a broad 
spectrum of Alfven waves, especially when the Alfvenic turbulence with phase coherent 
wave modes is given. Such particles may escape from the region of interaction with 
the Alfven waves and can contribute to fast particle population in astrophysical and 
space plasmas. 
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1 Introduction 



Plasma heating and acceleration are hot topics in the fields of nuclear fusion, plasma physics 
and astrophysics for a long time. A great number of heating mechanisms and theories 
have been proposed. In a collisionless space plasma environment, wave-particle interactions 
play significant roles of plasma heating and acceleration [l]-[6]. Even when collisions are 
included, plasmas can still be heated by the interplay between waves and particles to a 
certain degree |7-9|. Among a suite of electromagnetic waves, Alfven waves are generally 
thought to be the major contributor to the ion heating and acceleration since they exist 
pervasively in the solar wind and interplanetary space [lO] - [T3| . 

In thermodynamic sense, heating of ions by electromagnetic waves will lead to the dissi- 
pation of wave fields, thereby being an irreversible process. Several recent studies however 
show that ions can be heated by turbulent Alfven waves in low-beta plasmas even when 
no dissipation of wave fields occurrs a process that is called "pseudoheating" or 

nonresonant wave-particle interaction. Various methods have been used to validate this 
heating mechanism such as the test-particle approach with analytic solutions fl4|,[l7|, and 



quasi-linear theories 15,16 . In fact, the name "pseudoheating" may not be the best term to 



describe this "heating" process due to the kinetic effects of wave spectra [15, 19 . However, 
we still use this appellation for consistency in this paper. The pseudoheating is caused by 
the wave forces (or their spectra) that result in a deformation of the distribution function 
with respect to its initial Maxwellian shape. The mean-square velocity fluctuation due to 
the wave activity leads to an effective broadening of the distribution function similar to real 
heating and thus could mimic the genuine heating process [22]. It is important to point 
out that the heating shown by Ref. |14| contains both pseudoheating and genuine heating, 
which indicates that even in the particles' mean-velocity frame, the random kinetic energy 
of particles still increases via wave-particle interactions. The real heating (or the irreversible 
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dissipation of the wave fields) is caused by the initial pitch-angle scattering of newly created 
ions 1 20 1, indicating that there exists the dissipation of wave fields. In contrast, the first 
adiabatic invariant (magnetic moment /i = w±/B) remains constant during the pseudoheat- 



ing and thus it is an reversible process. The real heating in Refs. uA ) 20 is caused by the 
violation of the first adiabatic invariant due to the abrupt spatial change of the magnetic 
field. It is noteworthy that there is a factor of two difference between the temperature ex- 
pressions in Ref. [II] and Ref. [15], resulting from the fact that quasilinear theory cannot 
solve the problem of the wave damping appropriately. In addition, Wang et al. pointed out 
that pitch-angle scattering plays a key role in the pseudoheating process [l7| . 

In this paper, the heating process we focus on is restricted to the pseudoheating and thus 
the real heating is excluded. We demonstrate that the low- frequency Alfven wave propa- 
gating along the background magnetic field B = B i z can u heat" ions. The pseudoheating 
can be explained either by E x B drift (in the electric field of the Alfven wave and the 
ambient magnetic field B ) proposed in this paper or the pitch- angle scattering previously 
investigated [17]. Our analytic results, as will be shown below, are identical as those derived 
from quasilinear theory 15 . As shown in previous work [15, 17 : 



E^ = To f 1 + l^ 2 

n p °\Pp B 



T p± ~ T + — = T 1 + (1.1) 



where Wb = £>j^/2/x and n p are the wave magnetic field energy density and proton number 
density, respectively. (3 p denotes proton /3 that equals to (v) 2 jv\ ((v): the thermal speed). 
B^/Bq represents the ratio of wave-field energy density to that of the ambient field and 
T = m p (v) 2 /2fcg is the initial proton temperature. Here T p represents the " apparent tem- 
perature" that should be distinguished from the temperature associated with a real heating 
process [17] . 

The structure of the remainder of this paper is as follows: We derive and discuss the 
analytic results of pseudoheating based on E x B drift in Section 2. In Section 3, test particle 
simulation results are presented and discussed based on the comparison of a monochromatic 



Alfven wave and a wave spectrum. We also briefly discuss the recent observations of large 
amplitude magnetic perturbations associated with Alfven waves and their correlation to the 
pseudoheating. In the last section, conclusions are summarized. 



2 Analytic Theory of Pseudoheating 

Without loss of generality, left-hand circular polarized Alfven waves are considered in this 
paper. The wave magnetic field vector can thus be expressed as 

B w = ^2B k (cos(f) k i x -sin0 fe i y ), (2.2) 

k 

From the Faraday's law and the dispersion relation (Ampere's law) u = kv A) the electric 
field vector can be written as 

E w = -v A i z x B w (2.3) 

where \ X) i y and i z are unit directional vectors, (\) k = k(v A t — z) + tfk denotes the wave 
phase, (fk is the random phase for mode k and v A = B / \j\i^n p m p represents the Alfven 



speed. According to the linear approximation, we use Eq.(2.3) and the first order term of 
the generalized Ohmic law = — x B to derive the E x B drift velocity, v#, that can 
be expressed as follows: 



E w x B -v A i z x B w x B 



v £ = v ± = — 2 = — 2 (2.4) 

The derivation above is consistent with the traditional procedure to derive the classical 
Alfven wave solution [22] in the ideal magnetohydrodynamic (MHD) system. It indicates 
that the randomized proton motion is actually parasitic to wave fields due to the fact that the 
drift velocity is expected to disappear if the waves subside. Given the assumption that the 
characteristic spatial scale of the system is much larger than typical Alfven wavelength (lij , 
the energy density associated with the wave magnetic field can be written as: 

1 i d2 i d2 d2 d2 

t t t 1 // / \\2\ 1 2 D W 1 D D W D W /o x\ 

W B = -n p m p {(v E - {v E }) ) = -n p m p v A —^ = -n p m p — ^ = — (2.5) 



Bl 2 ^n p m p B$ 2/x 



where the bracket (•) denotes an average over all particles. Given the large characteristic 
spatial scale of the system, /, as described above, the following approximation is valid: 

ri 

l-HX> / J Q 



(v E ) = lim y / v E dz = (2.6) 



The temperature expression shown below is the same as Eq.(l.l) 



V = T + \n p m p ((v E - {v E )f) = T + || = T (l + l||) (2.7) 



The results shown in Eqs.(2.5)&(2.7) are consistent with the MHD theory due to the fact 
that the energy density associated with the Alfven wave magnetic field, £>j^/2// , equals 
to the ion (fluid) kinetic energy density, B^/2/hq = BqU^/ '(2fiov\) = p$u\/2, where U\ is 



the perturbed ion (fluid) velocity. The consistency between Eqs.(2.5)&(2.7) and the MHD 
theory I2TJ clearly indicates that the analytic theory described here incorporates the local 
equilibrium velocity distribution of the ions. In the following section, we adopt the test 
particle approach to simulate pesudoheating and the results will be presented and discussed 
in detail. 

3 Test Particle Simulations of Pseudoheating 

We start with a linearly polarized Alfven wave with wave magnetic field vector = 
J2k Bk cos (frkiy in order to show that the drift velocity is in E x B direction. This is the 
basis of further understanding. Then we conduct two case studies. In case one, we consider 
a monochromatic dispersionless Alfven wave with frequency oj=0.05Q p . In case two, we 
test a spectrum of Alfven waves with random phase and the frequencies of the wave 
modes can be calculated as follows: Ui = uj\ + (i — 1)Auj (i=l,2,..., A^; A=41), where 
Auj = (ujn — (jO\)/(N — 1); cji=0.01^ p and 6^=0.05^. The amplitude of each wave mode 
(only one wave mode in case one) is considered to be equal but changes gradually with time 
such that B^ = B\ = e{t)Bl, where 
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Figure 1: Wave field strength B^/Bq (solid line) and perpendicular apparent temperature 
T p j_ (dashed line) vs time; v p =0.07va- 



e(t) = 
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e e" 



-(t-tx) 2 /r 2 



if t<t 1: 

if h<t<t 2 , 
if t > t 2 . 



Similar to previous work [17|, we set ti=500fi~ 1 , i2=1000fi~ 1 , r=200^~ 1 , and 6o=0.05, 
where Q p is the proton gyrofrequency. The initial velocities of test particles are randomly 
distributed and possess a Maxwellian distribution with thermal speed v p =0.07va (the situ- 
ations with v p =0.01va, 0.03^, and 0.15va are also discussed later). The numerical scheme 
is similar to what described in Ref. (7, 14 . It is noteworthy that the ion thermal speed v p 

I 1/2 

in this paper and in Refs. 7 14 17 20 is denned as v p = (/c^T/m) 1 23 while the general 



thermal speed (v) = (2/c^T '/m) 1 ^ 2 \ causing v p = (v) 2 /2. Thus it leads to a factor of two 
difference for proton f5 p . The basic reason causing this difference is the different definition 



of temperature, i.e., T = m (v) 2 /2/cg or T = mv^/kB- The one-dimensional Maxwellian 
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Figure 2: Velocity scatter plots of test particles in the v x — v y space for a linearly polarized 
Alfven wave = ^ k Bj, cos at ^ p t=0 (left) and fi p t=600 (right); v p =0.07va- 



velocity distribution based on (v) and v p can be expressed as follows 23 



n f v? \ n ( v? 



f " ( w > = exp r J = (2^p exp <3 - 8) 

Different definitions of thermal speed or temperature, however, do not affect the final results 
since self-consistent definition is maintained throughout the previous studies (7([g[l5|[l7(|20|. 

Figjl] shows the dependence of apparent temperature (T p ±) on time-dependent wave 
field strength B^/Bq under a spectrum of Alfven waves. The result of a monochromatic 
dispersionless Alfven wave is almost the same as that under a spectrum of Alfven waves, 
which is in consistency with the analytic result shown above; the temperature expression 
Eq.(2.7) is independent of the number of wave modes, N. In FigjTJ apparent temperature 



versus time step shows the same tendency as the wave field strength, indicating that the 
proton temperature returns to its original value when the waves subside, in accordance with 
the work of Wang et al. [T7|[20]. It implies once again that the pseudoheating process is 
parasitic to the waves, as indicated by the aforementioned E x B drift. It is very important 



that the consistency between Eqs.(2.5)&(2.7) and the numerical solution shown in Fig{T] 
indirectly ensures the self-consistency of the test particle simulation debated in the Refs. 
[24p5] . Furthermore, this process is analogous to the ion motion in a magnetic mirror, where 
the magnetic moment is invariant. This physical picture helps us to better understand the 
reversibility in pseudoheating. 

In order to show that the pesudoheating is caused by E x B drift, we first illustrate the 
velocity scatter plots of test particles in the v x — v y space for a linearly polarized Alfven 



wave. According to Eq.(2.3), if is in the y direction, is in the x direction; therefore, 
E^ x B is in the y direction (ignore the negative sign here). As indicated in Figj2j the drift 
velocity is in the y direction, in agreement with the analytic results. To be consistent with 



the previous work 14-21 , we will focus on the circularly polarized condition in the following 
paragraphs. The further discussions are based on the comparison of a monochromatic Alfven 
wave and a wave spectrum. 

FigJS] presents the velocity scatter plots of ions for a circularly polarized Alfven wave 
with different initial thermal speeds: v p =0.01va, 0.03^, 0.07va and 0.15va> Inspection 
of FigjS] reveals that when the thermal speed is quite small with respect to the Alfven 
speed va (i.e., v p =0.01va and 0.03^), test particles form a ring distribution in the v x — v y 
velocity space. Although the ring distribution eventually can be filled with test particles 
with different velocities when v p becomes large, the heating efficiency becomes relatively low 



[refer to Eq.(2.7)]. Besides the velocity scatter plot, particle distribution in the phase space 
is also essential for our understanding of the pesudoheating. Fig|4] shows the scatter plots 
of protons between 1000 and 2400^ A ^- 1 at different times Q p t=0, 200, 300, 600, 1200 and 
1500 for a monochromatic dispersionless Alfven wave. The results agree with that shown 
in Figjl| the stronger the wave fields are, the more obvious the velocity fluctuations are. 
A common feature that stands out in Figj4] is that the particle motion is periodic under a 
monochromatic Alfven wave. It indicates that the kinetic behavior of test particles under a 
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v =0.157, 
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Figure 3: Velocity scatter plots of test particles in the v x — v y space for a circularly polarized 
Alfven wave at f2 p t=600. (a) v p =0.01va (b) v p =0.03va, (c) v p =0.07va, and (d) v p =0.15va- 
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Figure 4: Scatter plots of protons between 1000 and 2400^ 1 at different times, fi p t=0, 
200, 300, 600, 1200 and 1500 in the v x -z phase space for a circularly polarized Alfven wave; 
v p =0.07va> 

monochromatic dispersionless wave versus a wave spectrum, as will be shown below, is fairly 
different. 

Compared with the ion behavior under a monochromatic wave, the ion motion tends 
to become random under a wave spectrum as indicated in Figsj5] & [6j Figj5] presents the 
velocity scatter plots of test particles in the v x — v y space for a spectrum of circularly 
polarized Alfven waves with random phases ipk- When the thermal speed is two orders of 
magnitudes smaller than the Alfven speed va (i.e., v p =0.01va and 0.03^), test particles 
cannot fully fill the circle in the v x — v y velocity space. However, with the increase of the 
initial thermal speed v p , the circle in the velocity space is fully filled with test particles with 
different velocities. Figsj5](e)&(f) show the velocity scatter plots when adopting 1001 wave 
modes. Compared with Figsj5](a)&(b), Figsj5](e)&(f) show that protons tend to fully fill 
the circle in the velocity space, indicating that the ion velocity distribution tends to be a 
full-filled circle in the velocity space when the number of wave modes, AT, is large enough, 
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Figure 5: Velocity scatter plots of test particles in the v x — v y space for a spectrum of 
circularly polarized Alfven waves with random phases cpk at ^ p t=600. (a) v p =0.01va, (b) 
v p =0.03va, (c) v p =0.07va, (d) v p =0.15va, (e) v p =0.01va, and (f) v p =0.03va', ( a )-(d): 41 
wave modes, (e)-(f): 1001 wave modes. n 



regardless of the relatively small thermal speed. It can be observed based on the comparison 
between Figj3] and Figj5] as well. The phase space proton scatter diagrams by adopting a 
wave spectrum with wave modes 7V=41 are shown in FigjH] . The results, however, reveal 
significantly different particle distributions in the phase space compared with those shown 
m Figg} The main conclusion drawn from FigjH] is that the velocity fluctuations caused 
by wave activity is quasirandom, and thereby could mimic the real heating (also refer to 
Figj8]). However, as indicated in both Fig j2] and Figj6| the pseudoheating caused by these 
wave activities is reversible, indicating no dissipation of wave fields, and therefore does not 
represent real heating in thermodynamic sense. It is also interesting to investigate the ion 
behavior under a spectrum of circularly polarized Alfven waves with same initial phases (fk- 
It is noteworthy that the wavelength and wave frequency among different wave modes are 
still different. Fig{7] illustrates the scatter plots of the test particles in the v x — z and v x — v y 
space. There is a pulse-like structure in the protons' phase space distribution due to the 
coherence of wave modes. The proton distribution in the v x — v y space is primarily consist 
of two parts: the broadening core Maxwellian distribution and the outer ring structure. The 
broadening is caused by the wave forces (or their spectra) while the accelerated particles 
in the outer ring result from the Alfvenic turbulence with phase coherent wave modes as 
indicated in Ref. 126), where the acceleration of charged particles by large amplitude MHD 
waves was studied. In contrast, if Alfvenic turbulence without any envelope modulation (27] 
is given, the "acceleration" may not be observed. These high energy particles in the outer 
ring may escape from the region of interaction with the Alfven waves and can contribute to 
the fast particle population in astrophysical and space plasmas. 

In FigjSJ we present the normalized velocity distribution function at different time Q p t=0 
(solid line) and 600 (dashed line). The broadening of the distribution functions shown 
in Figsj8](a)-(c) are based on the statistics of all the particles with different spatial co- 
ordinates z while FigMd) is based on the local statistics of protons in the spatial range 



12 



(a) 



■(b) 




1000 1200 1400 



600 1800 2000 2200 2400 
ZD /V 



1000 1200 1400 



600 1800 2000 2200 2400 
ZD /V 






1000 1200 1400 1600 



1800 2000 
ZD /V 



2200 2400 



Figure 6: Scatter plots of protons between 1000 and 240(huf2~ 1 at different times, f2 p t=0, 
200, 300, 600, 1200 and 1500 in the v x -z phase space for an Alfven wave spectrum with 
random phases cp^ v p =0.07va- 




Figure 7: Proton scatter plots of a spectrum of circularly polarized Alfven waves with same 
initial phases tpk at fi p t=600; v p =0.07va- 
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Figure 8: The normalized velocity distribution functions plotted against v x at Q p t=0 (solid 
line) and 600 (dashed line), (a) A spectrum of circularly polarized Alfven waves with ran- 
dom phases (fk] (b) A circularly polarized monochromatic Alfven wave. v p =0.07va', (c) A 
spectrum of circularly polarized Alfven waves with same initial phases (fk- (d) The local 
normalized velocity distribution function based on the statistics of protons in the spatial 
range 1490 < zQ p /va < 1510 under a spectrum of circularly polarized Alfven waves with 
random phases ipk- The black dots denote the shifted local normalized velocity distribution 
function in the laboratory frame while the white dots represent the local normalized veloc- 
ity distribution function in the particles' mean-velocity frame. The local statistical mean 
velocity equals to ^(1490 < zQ p /va < 1510) = — 0.176^ • 
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1490 < zQ^p/ v a < 1510, thus being a local normalized velocity distribution function around 
zQ p /va = 1500 as indicated by the line in Figj6](d). The velocity spreading in Figsj8](a)-(c) 
is caused by averaging over the wave effects. Inspection of Figsj8](a)-(b) reveals that the 
Maxwellian distribution is more likely to maintain under a wave spectrum compared with a 
monochromatic Alfven wave due to the fact that the ion motion under a wave spectrum is 
quasirandom, which is in good agreement with the analytic derivation by Wu and Yoon using 
the quasi-linear theory PL5]. It also points towards different kinetic behaviors of particles 
under a wave spectrum and a monochromatic Alfven wave. It is well known that enhanced 
Alfvenic turbulences exist pervasively in the solar wind and interplanetary space, while in 
astrophysical observation measurements of temperature are based on spectroscopic data col- 
lected from the source region of interest. Therefore it is very difficult to distinguish the real 
heating from the pseudoheating due to the restriction of spatial resolution of the instruments 
and the presence of wave forces (or their spectra). The suprathermal tail and small bump 
shown in Figj8](c), where the v x axis scale is different from Figsj8](a)-(b), corresponds to the 
accelerated (suprathermal) ions due to the wave modes coherence. The bump-on-tail struc- 
ture may excite plasma instabilities, however, the detailed discussion of these instabilities is 
beyond the scope of this paper. It also needs to point out that here we only consider the 
pseudoheating and investigate the normalized velocity distribution function corresponding 
to this process while the real heating via nonresonant interaction with Alfven waves always 
coexists [l4j[25]. The velocity distribution function, therefore, is supposed to be broader 
than that shown in FigjS] and it is possible that the broadening effects smooth the velocity 
distribution function and eliminate the bump on the tail. In Figj8](d), we investigate the 
local normalized velocity distribution function based on the statistics of protons in the spa- 
tial range 1490 < zQ p /va < 1510. If the spatial range is too small, there will be insufficient 
particles to be counted for statistics. The local statistical mean velocity equals to v x =- 
0.176^ which agrees well with the result shown in Figj6](d). The local normalized velocity 
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distribution function in the particles' mean-velocity frame at Q p t = 600 (as indicated by the 



white dots) is in good agreement with the initial Maxwellian distribution based on Eq.(3.8). 
The slight difference between these two distribution functions results from the fact that the 
local statistics are based on a small spatial range (1490 < zQ p /va < 1510) and therefore it 
can be potentially affected by the effects of wave spectra. This leads to a slightly broader 
local velocity distribution. The result shown in Figj8](d) indicates that the analytic theory 
in Sec. 2 is on the basis of the local equilibrium velocity distribution of the ions. 



The pseudoheating becomes important when the magnetic perturbation (SB) is rela- 
tively large compared with the background magnetic field (Bo). Figj9] illustrates the spatial 
distribution of the Alfven wave magnetic fields with different wave modes, TV, and the cor- 
responding scatter plot in the v x — z phase space. The solid curve for the magnetic field of 
a monochromatic Alfven wave is simply a sinusoidal wave. The dashed and dotted curves 
represent a spectrum of Alfven waves with 41 wave modes that has random and same phase 
(fki respectively. They can mimic the Alfvenic turbulence due to the quasirandom fluctu- 
ations of the wave fields. The amplitude of the magnetic perturbation can be quite large 
locally, thus may produce some energetic particles that can escape from the constraint of 
the Alfven waves and contribute to fast particle population in astrophysical and space plas- 
mas. As indicated in Figj9j the larger the magnetic perturbations are, the more effective 



the distribution functions are broadened. Refs. |32, 33 suggest that small scale reconnection 
events occur during the solar flares, which can provide large magnitude spike-like magnetic 
field fluctuations. Additionally, the cluster observation of surface waves in the ion jets from 
magnetotail reconnection also shows that SB/\B\ can be as high as 0.5 and occasionally 
even higher |34|. Furthermore, the amplitude of the magnetic perturbation associated with 
Alfven waves can be even larger than the ambient magnetic field in certain cases observed 
by the Wind satellite [35] . Figsj9](b)&(c) also show that the distribution of test particles in 
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Figure 9: (Color online) (a) The spatial distribution of Alfven wave magnetic fields at 
^ p t=600. (b)-(c) The Alfven wave magnetic field spatial distribution and the corresponding 
scatter plot in the v x — z phase space at ^ p t=600; (b) the case with random phases (fk and 
(c) the case with same phases cp^. 
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the v x -z phase space and the wave field perturbations are in antiphase (tt phase difference) 
due to the fact that the proton motion is parasitic to Alfven waves, indicating an exchange 
of energy between the particles' kinetic energy and the magnetic energy. Therefore Figj9](b) 
indicates that the psuedoheating is a consequence of equilibrium MHD system. 

It is noteworthy that the parallel electric field and wave damping are not included in the 
test particle simulations of our current work. In the self-consistent simulations such as hybrid 
simulations (i.e., kinetic description of ions, fluid electrons), the phase-correlation among 
the Fourier modes could affect, for instance, that the ponderomotive force resulting from 



the envelope-modulated Alfven waves heat ions through the nonlinear Landau damping 28 



and the parallel heating of ions due to the nonlinear Landau damping [29]. In Ref. [28] , 
the ponderomotive force and beat interaction are identified as the most important nonlinear 
effects in proton heating by nonlinear field-aligned Alfven waves in the solar coronal holes. 
Interestingly, they found that the nonlinearity is particularly strong when the wave spectrum 
consists of counterpropagating modes of equal intensity, even if the intensity is relatively low. 
Moreover, from the hybrid approach, dissipation processes of the Alfvenic turbulence with 
the broadband wave number spectrum can be different from those of monochromatic Alfven 
waves, since the former is associated with the density fluctuations, \5n/n\, and the resultant 



spatial modulation of \B\ due to compressive effects of ponderomotive forces [30,31 : right- 
hand polarized Alfvenic turbulence with such an envelope modulation can be dissipated due 
to the nonlinear Landau damping, while left-hand polarized Alfvenic turbulence with the 
broadband spectrum is preferentially dissipated by the modulational instability. On the 
other hand, low-frequency, monochromatic Alfven waves are relatively stable to the linear 
collisionless dissipations such as Landau damping and cyclotron damping, nonlinear wave- 
wave interactions parametric instabilities are important for the dissipation of these waves. 
The decay instability is dominant for dissipation of the right-hand polarized finite-amplitude 
monochromatic Alfven waves in low /3 plasmas, while the left-hand polarized waves can 
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also be dissipated via the modulational instability 29 . These wave damping mechanisms 
are closely related to the present work; therefore, in order to study the characteristics of 
the "acceleration" and "heating" process in more detail, comprehensive and self-consistent 
studies are necessary in the future. 

4 CONCLUSION 

Ion pseudoheating by low-frequency Alfven waves is investigated based on the comparison 
between a monochromatic Alfven wave and a wave spectrum. Both analytic and simulation 
results show that E x B drift plays a principal role in this process and the proton motion is 
parasitic to Alfven waves. It indicates the psuedoheating is a consequence of the equilibrium 
in the MHD system. Our results are in good agreement with the previous studies based on 
pitch-angle scattering. More importantly, it provides a simple understanding of the reversible 
property of this process from E x B drift point of view; if wave magnetic and electric fields 
disappear, there will be no drift velocity ve and therefore no pseudoheating. We showed 
that the wave spectra contribute to the broadening of the Maxwellian distribution during 
the pseudoheating, and it is therefore difficult to exclude the apparent temperature T p from 
observations due to the low spatial resolution of the instruments. It is of particular interests 
to note that the Maxwellian shape is more likely to maintain during the pseudoheating under 
a wave spectrum compared with a monochromatic Alfven wave. We, therefore can conclude 
that the kinetic behavior of ions under a monochromatic wave and a wave spectrum is totally 
different. Moreover, we illustrated that E x B drift can produce energetic particles under a 
spectrum of Alfven waves, which may contribute to fast particle population in astrophysical 
and space plasmas. 
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